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ABSTRACT 

The  world  of  vehicle  design  is  a  fast-paced  iterative  environment  that  demands  efficiency  in  the  simulation  of 
suspension  loads.  Toward  that  end,  a  computationally  efficient,  linear,  planar,  quasi-static  tire  model  is  developed  in 
this  workthat  accurately  predicts  a  tire ’s  lower-frequency,  reasonably  large  amplitude,  nonlinear  stiffness  relationship. 
Hamilton’s  principle  is  used  to  derive  the  axisymmetric  and  circumferentially  isotropic  stiffness  equation  which  is 
discretized  into  segments  to  create  the  tire  stiffness  model.  The  model  is  parameterized  by  a  single  stiffness  parameter 
and  two  shape  parameters  such  that  the  tire’s  deformed  shape  is  completely  independent  of  the  overall  tire  stiffiiess 
and  the  forces  acting  on  the  tire.  Constraint  modes  are  developed  that  capture  the  tire  enveloping  and  bridging 
properties  via  component  mode  synthesis  originated  by  Hurty  and  Gladwell  and  the  Guyan  static  reduction  method 
decouples  active  constraints  fromthe  stiffiiess  matrixformulation.  A  recursive  method  is  developed  to  deduce  the  set 
ofactive  constraints  at  the  tire -road  interface.  The  model  captures  the  nonlinear  stiffiiess  of  a  real  tire  by  enforcing 
the  unidirectional  geometric  boundary  conditions  during  the  recursive  method.  The  model  parameters  are  identified 
via  two  quasi-static  experiments:  a  flat -plate  and  a  cleat  test.  The  simulated  vertical  loads  are  within  4%  of  the 
experimental  load  throughout  a  reasonably  large  range  of  travel  for  the  flat -plate  test  and  within  7%  of  the  experimental 
loads  for  the  cleat  test.  The  simulation  produces  nonlinear  stiffiiess  when  simulating  a  flat  plate  test  and  a 
discontinuous  stiffiiess  when  simulating  a  cleat  test.  This  model  strikes  a  balance  between  simple  tire  models  that 
lack  the  fidelity  to  make  accurate  chassis  load  predictions  and  computationally  intensive  models  that  cannot  provide 
timely  predictions.  It  is  expected  that  this  tire  model  is  used  as  an  integral  part  of  computationally  efficient  and 
accurate  vehicle  simulations  that  are  critical  throughout  the  fast-paced  iterative  design  process. 
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NOMENCLATURE 

N  Totalnumberoftire  segments 
u  Set  of  radial  deflections 
R  Tire  radius 

Fn  Radial  force  acting  on  nth  tire  segment 
Fa  Set  of  active  radial  tire  force 
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F  Vertical  spindle  force 
/  Linearly  distributed  load  density 
k  Radial  stiffness  density 
E  Elastic  modulus 
/  Second  moment  of  inertia 

G  Shear  modulus 
A  Cro  s  s  -s  ectional  area 

Ub  Distributed  bending  potential  energy 
Us  Distributed  shear  potential  energy 
Ur  Distributed  radial  elastic  potential  energy 
Uw  Distributedextemalwork 
yb  Distributed  bending  stiffness 
ys  Distributed  shear  s tiffh ess 
yr  Distributed  radial  stiffness 
k  o  Model  s  tiffh  ess  parameter 
a1a2  Model  shapeparameters 


INTRODUCTION 

With  the  advent  of  autonomous  vehicles  and  the  simultaneous  pressure  to  improve  energy  efficiency,  there  aie  ever 
increas  ing  demands  on  vehicle  engineers.  In  particular,  chassis  engineers  are  faced  with  numerous  design  challenges 
(weight,  cost,  ride,  handling,  noise,  reliability,  packaging...)  while  interacting  with  several  influential  groups 
(scientific  labs ,  proving  ground,  suppliers ,  safety,  manufacturing,  powertrain,  body. . . )  that  result  in  a  highly  iterative 
design  process.  All  these  challenges  and  groups  must  be  satisfied,  to  the  best  ofthe  engineer’s  ability,  within  a  veiy 
short  timeframe.  This  fast-paced  iterative  design  environment  demands  efficiency.  Yet  some  fundamental  physical 
truths  remain  within  this  cacophony  of  ever-changing  requirements,  influential  groups,  and  deadlines.  For  example, 
regardless  of  how  the  power  is  developed,  or  how  the  wheels  are  steered,  all  the  vehicle  forces  must  ultimately  be 
reacted  through  the  vehicle  suspension  andtires . 

The  goal  then,  is  to  create  a  highly  efficient  model  that  captures  the  fundamental  physical  properties  of  tires  from 
which  simulation  results  are  used  to  make  informed  design  decisions.  High  fidelity,  but  computationally  intensive 
models  that  provide  very  accurate  information  after  the  parts  have  been  released  for  production  are  of  no  use. 
Simulation  results  arising  fromoverly  simplified  models  that  lack  the  fidelity  to  make  accurate  chassis  load  predictions 
cannot  be  trusted  and  have  little  value.  The  objective  of  this  work  is  to  develop  a  computationally  efficient,  planar 
tire  model  that  accurately  predicts  the  lower-frequency,  but  not  necessarily  low  amplitude,  tire  shape.  Higher- 
frequency  dynamics  such  as  those  thatwouldbe  required  for  noise  predictions  are  outside  the  scope  ofthis  work  The 
emphasis  is  intentionally  placed  on  the  shape  ofthe  deformed  tire  in  the  region  of  the  contact  patch,  which  must  be 
capable  of  representing  the  bridging  property  (over  narrow  cracks)  and  enveloping  property  (over  sharp  bumps).  The 
scope  includes  relatively  large  deformations,  approaching,  but  not  including,  the  point  of  rim  strike.  Although  the 
relationship  between  the  radial  deformation  and  radialload  is  developed  as  a  linear  quasi-static  model,  the  nonlinear 
stiffness  relationship  between  the  vertical  spindle  load  and  the  vertical  tire  deflection  that  exists  in  real  tires  is  captuied 
by  the  model  via  the  unidirectional  geometric  boundary  conditionat  the  tire -road  interface.  Six  physical  tire  properties 
(radius,  bending  stiffness. . .)  are  reduced  in  the  development  of  this  model  to  a  single  stiffness  parameter  and  two 
shape  parameters.  Experimental  data  are  used  to  validate  the  approach  using  two  quasi-static  tests,  a  flat  plate  and  a 
cleat  test,  with  relatively  large  deformations. 

The  remainder  of  this  work  is  developed  as  follows.  Some  background  on  finite  element  and  finite  difference  tie 
modeling  and  the  use  of  constraint  modes  in  tire  modeling  is  briefly  reviewed.  The  mathematical  derivation  of  the 
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parameterization  of  the  stiffness  matrix  is  developed  followed  by  the  development  of  the  constraint  mode  model  for 
this  work.  A  recursive  process  to  determine  the  proper  tire  shape  is  developed  followed  by  validation  with  quasi¬ 
static  experimental  results.  The  resulting  mo  del  is  capable  ofbridging  and  enveloping  irregularities  in  the  roadsurtace, 
generating  the  correct  tire  shape  over  irregular  surfaces,  and  capturing  the  nonlinear  stiffness  relationship  produced 
by  physical  tires. 


BACKGROUND 

One  long-standing  hurdle  to  accurately  predicting  responses  (forces,  torques,  and  accelerations)  at  the  tire  spindle  is 
the  development  of  high  fidelity  tire  models.  A  considerable  amount  of  tire  modeling  research  has  been  conducted 
and  a  comprehensive  review  of  the  current  stateofthe  art  is  presented  by  Willumeit  and  Boehm  [1]  and  more  recently 
by  Umsrithong  [2].  Although  a  complete  review  of  recent  developments  is  not  presented  herein,  the  results  horn 
several  of  these  works  have  direct  relevance  to  the  objectives  of  this  study. 

Researchers  using  Finite  Element  (FE)  models  have  demonstrated  that  spindle  responses  can  be  predicted  accurately, 
but  may  be  computationally  intensive  and  require  a  great  number  of  physical  parameters  that  are  difficult  or  costly  to 
obtain.  Mousseau  andHulbert  [3]  and  Darnell,  Hulbert,  and  Mousseau  [4]  developed  a  3-D  tire  modelusing  an  FE 
modelfor  the  tread  anda  lookup  table  for  the  sidewall  membrane  elements.  Hanley,  Crolla,  and  Hauke  [5],useda  1- 
D  FE  model  for  tires  traversing  large  obstacles;  a  more  computationally  efficient  means  to  simulate  relatively  laige 
deflections  is  pursued  in  the  present  study.  Burke  and  Olatunbosun  [6]  used  an  FE  model  and  a  contact  patch 
description  to  predict  vertical  force  and  displacement.  Presently  a  more  formal  approach  to  defining  the  contact  patch 
is  developed  through  the  definition  of  the  active  constraints  between  the  tire  and  road  surfaces.  The  current  state  of 
FE  modeling  work  is  presented  by  Ragheb,  El-Gindy,  and  Kishawy  [7]  wherein  Mooney-Rivlin  material  properties 
were  used  to  create  detailed  FE  models  of  the  tread  and  tire.  Their  quasi-static  flat-plate  test  results  are  used  as  a 
baseline  to  validate  the  capability  ofthe  work  developed  in  the  present  study. 

Other  researchers  have  developedring  models  oftires  using  finite  differences.  Van  Oosten  and  Jansen  [8]  useda  rigid 
ring  to  calculate  belt  vibrations  while  the  deformed  shape  ofthe  tire  (including  bridging  and  enveloping)  was 
accurately  predicted  by  Zegelaar  and  Pacejka  by  using  a  flexible  ring  to  simulate  the  quasi-static  response  of  a  tite 
rolling  over  an  uneven  surface  [9].  Although  the  simulation  results  agreed  with  the  experimental  data,  the  model  is 
computationally  intensive  andrequires  comprehensive  physical  testing  to  define  the  parameters.  Asimilar  approach 
was  developedby  Loo,  who  modeled  the  tire  by  using  a  more  comp  lex  flexible  circular  ring  under  tension  with  a  nest 
oflinearspringsanddampersarrangedradially  [10].  However  in  Loo’s  study  thering  tension  and  foundation  stiffiless 
requires  additional  contact  patch  measurements  and  the  parameters  were  experimentally  identified  using  multiple 
single-point  load  tests.  Badalamenti  et  al.  showed  that  a  radial  spring  tire  model,  in  which  the  radial  spring  element 
deflection  depends  on  the  adjacent  element  deflections,  could  accurately  describe  the  tire  enveloping  behavior  in  an 
efficient  model  [1 1].  Presently,  a  planar  ring  model  with  an  elastic  foundation  that  includes  shear  and  bending  stiffiiess 
in  ring  is  developed  to  accurately  and  efficiently  describe  the  tire  enveloping  behavior. 

It  is  assumedin  this  workthat  the  low-frequency  deformation  ofthe  tire  shape  can  be  adequately  defined  by  the  quasi¬ 
static  constraintmodes.  This  as  sumption  is  given  credence  by  Gillespie’s  studyin  which  a  radial  spring  mo  del  is  used 
to  simulate  the  stiffiiess  variation  circumferentially  and  the  magnitude  of  the  radial  force  variation  is  found  to  be 
relatively  independent  of  speed  [12].  Similarly,  Takayama  et  al.  developed  a  model  to  predict  the  transient  icsponse 
of  a  tire  encountering  a  cleat,  where  the  belt  and  tread  region  is  modeledby  a  rigid  ring  and  deflections  fromthe  cleat 
are  absorbed  by  a  linear  and  planar  spring  attachedto  the  rigid  ring  [13]. 

Ferris  suggests  that  a  static  constraint  mode  tire  model  can  be  used  to  capture  the  tire  enveloping  and  bridging 
properties  with  an  axisymmetric  and  circumferentially  isotropic  model  [14].  This  reduced  representation  of  the  tite 
model  is  founded  on  the  seminal  component  mode  synthesis  work  originated  by  Hurty  [15]  and  Gladwell  [16]. 
Presently,  in  Equations  (15)-(18),  the  order  ofthe  stiffiiess  matrix  is  reduced  via  a  Guyan  static  reduction  [17].  This 
constraint  mode  tire  model  was  expandedby  Ma  to  include  shear  effects  in  the  derivation  and  to  reduce  the  requited 
number  of  parameters  fromfive  dependent  parameters  to  two  parameters  that  define  the  constraintmodes  (theshape) 
and  an  independent  parameter  that  defines  the  overall  tire  stiffiiess  [18].  Presently,  these  concepts  are  combined  to 
develop  a  complete  formulation  ofthe  planar  constraint  mode  tire  model. 
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MODEL  DEVELOPMENT 


Consider  the  tire  to  be  a  ring  that  is  inextensible  in  total  circumference,  but  radially  flexible,  and  supported  by  an 
elastic  foundationparameterized  by  the  radial  stiffh  ess,  k,  as  showninFigure  la.  The  assumed  axis  ymme  trie,  plan  ai; 
and  circumferentially  isotropic  tire  is  divided  into  TV  segments.  Each  segment  is  modeled  as  an  Euler  elastic  beam  of 
length  2  71R/N  and  notional  cross-sectional  area,  A,  shown  as  a  shadedregion  in  Figure  lb  (the  model  itself  is  planar). 
The  vector  {u}  defines  the  radial  deflection  of  all  tire  segments  and  the  quasi-static  constraint  modes  are  devebped 
from  the  stiffness  relationship  given  in  Equation  (1). 


Figure  1.  Schematics  of  (a)  flexible  ring  with  radial  stiffness  and  (b)  single  tire  segment 


[K]{u}={F}  (1) 

Note  that  this  derivation  includes  relatively  large  deformations,  approaching,  but  not  including,  thepoint  of  rims  trite 
(when  the  tire  sidewall  is  folded  such  that  the  upperportion  is  in  direct  contact  with  the  lowerportion  and  there  is  a 
discontinuous  and  dramatic  increase  in  apparent  stiffh  ess).  That  is,  the  radial  stiffh  ess  of  each  segment  is  assumed  to 
be  linear.  Next  the  stiffness  matrix[7T|  is  developed,  followedby  a  reduced  order  parameterization  of  the  matrixusing 
one  stiffness  andtwo  shape  parameters. 


Derivation  of  the  Stiffness  Matrix 


The  quasi-static  equation  resulting  fromthe  application  ofHamilton’s  principle  forthe  distribution  o  f  energy  and  woik 
throughout  the  tire  is  written  as  Equation  (2). 

ft2  (2) 

(8Ub-8Us-8Ur  + 8Uw)dt  =  0  Vtlft2 

The  virtual  changes  are  computed  as  Equations  (3)  -  (6),  where  spatial  derivatives  are  represented  by  the  commonly 
used  prime  notation  (e.g.,  first  spatial  derivative  ofradial  displacement  is  u ') 


8Ub  =  EIvT'Su  (3) 

8US  =  GAu"8u  (4) 

8Ur  =  kuSu  (5) 

8UW  =  f8u  (6) 


Continuity  is  enforced  through  the  boundary  conditions  on  the  radial  deformations  and  their  spatial  derivatives  as 
given  in  Equation  (7). 


x  ^\x+2nR>  U  \x  U  \x+2nR>  ^  lx  ^  I 


x  +  2nR> 


u'"\x  =  W"  | 


x  +  2nR 


The  stiffness  equation  for  each  segment  is  then  simplified  to  Equation  (8). 


(7) 


EIu""  +  GAu"  +  KU  =  f 


(8) 


To  reduce  the  required  notation,  the  sixphysical properties  (E,  I,  G,  A,  k,  and  R)  are  grouped  into  three  tire  stiffiiess 
parameters  (yb,  ysandyr)  as  defined  in  Table  1.  The  tire  stiffiiess  parameters  are  invariant  for  a  specific  tire  and 
positive  valued  for  all  tires .  Mostnotably  the  tire  stiffiiess  parameters  are  not  a  function  of  the  discretization  process. 
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Table  1.  Tire  stiffness  parameters  and  discretized  radial  force 


yb  =  EI(2nRy3 

Distributed  bending  stiffiiess 

Ys  =  GA(2nR)~1 

Distributed  shear  stiffiiess 

yr  =  k(2ttR) 

Distributed  radial  stiffiiess 

Fn  =  f-(2nR) 

Radial  force  on  nth  tire  segment 

The  stiffness  equation  for  each  segment  is  then  rewritten  as  Equation  (9). 

yb(2nR)3u""  +  ys(2nR)u"  +  yr(2nR)~1u  =  f  (9) 


Defining  Stiffness  and  Shape  Parameters 

Next  the  model  is  discretized  into  N  equal  segments  via  a  finite  difference  method  The  fourth  spatial  derivative  of 
the  radial  displacement  for  the//7  tire  segment  is  approximated  by  Equation  (10);  the  second  spatial  derivative  is  found 
in  a  similar  manner. 


u 


rrrr  ^ 

n  ~ 


^n-2  4un_1  +  6un  4'Un+1  +  'Un+ 2 


(10) 


A  second  set  of  parameters  is  introduced  to  simplify  the  subsequent  analysis,  beginning  with  discretizing  the 
continuous  stiffness  equation  given  in  Equation  (9).  The  discrete  stiffness  equation  for  the  nth  tire  segment  is  expressed 
as  Equation  (11),  where  the  radial  force  on  the  //tire  segment,  is  defined  in  Table  1. 

k  Qa2lLn_2  "h  ^Oalun-l  "h  ^O^T^n+1  "b  ^0^2 ^n+2  Zi 


The  model  stiffiiess  (ko)  and  shape  parameters  (ai  and  ai)  are  defined  in  terms  of  the  tire  stiffness  parameters  of  Table 
1.  For  convenience,  the  conversion  between  sets  ofparameters  is  defined  by  a  pair  of  transformation  matrices  given 
in  Equation  (12)  and  Equation  (13). 


1  1 
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0 

N3 
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7C 
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72 
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.IV3  0  0. 

u  5  A 

(13) 


The  usefulness  of  the  tire  stiffiiess  parameters,  {y),  is  that  they  are  convenient  to  describe  the  distributed  stiffiiess 
equation,  Equation(9),  and  they  are  notaffectedby  the  number  of  tire  segments  chosen  to  discretize  the  system.  The 
model  stiffiiess  and  shape  parameters  (ko,  ai ,  ai),  by  contrast,  are  convenient  for  the  discretized  stiffiiess  equation, 
Equation  (1 1),  with  which  the  reminder  of  this  work  is  concerned.  Consider  a  typical  modeling  exercise  in  which  the 
number  of  segments  is  to  be  increased.  Once  the  model  stiffiiess  and  shape  parameters  are  identified  for  a  model  with 
the  original  number  of  segments,  the  corresponding  invariant  tire  stiffiiess  parameters,  {y},  can  be  calculated  using 
Equation  (12).  If  a  new  number  of  segments  is  desired,  the  corresponding  values  for  the  new  model  stiffiiess  and 
shape  parameters  is  calculated  usingEquation  (13)  and  the  invarianttire  stiffiiess  parameters. 

Parameterized  Stiffness  Matrix 

The  symmetry  of  a  physical  tire,  enforced  by  Equation  (7),  results  in  a  circulant  Toeplitz  stiffiiess  matrix  [K],  as  shown 
in  Equation  (14),  that  comprises  two  parts:  the  single  parameter  k0  defines  the  model  stiffiiess  and  the  model  shape 
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matrix  [a]  (defined  by  the  two  model  shape  parameters  ai  and  ai)  that  defines  the  relative  deformations  of  the  tie 
when  acted  upon  by  geometric  constraints. 


DEVELOPMENT  OF  CONSTRAINT  MODES 

The  Guyan  static  reduction  method  [17]  is  used  to  reduce  the  order  of  the  stiffness  matrix  in  Equations  (15)-(18)by 
imposing  geometric  boundary  constraints  on  the  active  degrees  of  freedom.  The  degrees  of  freedom,  defined  as  the 
radial  deformations  of  the  tire  {w},  are  categorized  into  active  and  omitted  sets,  representedby  superscripts  ‘<2 ’and  ‘o’ 
respectively.  In  general,  the  geometric  boundary  constraints  mustbe  a  subset  of  the  active  degrees  of  freedom.  In  the 
tire  modeling  cas  e  with  which  this  work  is  concerned,  the  active  constraints  are  exclusively  compris  ed  of  those  degiees 
of  freedomin  contact  with  the  road  surface.  The  degrees  of  freedomare  reordered  and  written  as  vector  { u }  such  that 
the  active  degrees  of  freedomoccupy  the  first  positions  in  the  vector,  followed  by  the  omitted  degrees  of  freedom,  via 
an  orthogonal  sorting  matrix,  S,  as  shown  in  Equation  (15). 

[u}=[[Sa]  [5°]]{[J]j=  [S]{u],  {u}=[S]r{u}  (15) 

The  stiffness  matrix  and  the  force  vector  must  be  similarly  reordered  as  given  in  Equation  (16). 

[£]  =  [S]r[tf][S],  {f}=  [5]t{F]  (16) 

The  reordered  stiffness  equation  is  then  partitioned  according  to  the  active  and  omitted  degrees  offreedomas  given 
in  Equation  (17). 

»«=[[£ 

Given  a  set  of  active  constraints,  [ua  },  the  omitted  constraint  vector  {u°}  is  written  as  Equation  (18).  The  constraint 
mode  for  this  particular  set  of  active  constraints  is  written  as  Equation  (20). 

{1 u° }  =  -[^°°]'1[^oa]{ua}  (18) 

In  this  work,  the  stiffness  matrix  is  defined  in  terms  of  an  overall  stiffness  parameter  and  two  shape  parameters,  as 
developed  in  Equation  (14).  The  reordered  stiffness  submatrices  are  then  givenby  Equation  (19)  as 

[K°°]  =  [S°r  Vc0[a])[S°]  and  [R°a]  =  [5°]r  (fc0[a])[5“]  (19) 

The  omitted  constraint  vectorbecomes  solely  dependent  on  the  shape  parameters  and  theactive  constraint  vectoi;  as 
developed  in  Equation  (20). 

{u0}  =  -([Stf],,[a][S0])-1([S0]I,[a][Sffl]){fi“}  (20) 

Finally,  the  complete  displacement  vector  is  written  as  a  function  ofthe  active  constraint  vector  in  Equation  (21). 

{u}  =  [[S“]  -  [S°]  ([S°]T  [a]  [S°]  r1  ([5°] r [a]  [S“] )]{«“}  (21) 

Note  that  this  mode  shape  is  completely  independent  of  the  overall  tire  stiffness  and  the  forces  acting  on  the  tiie. 
Rather,  the  shape  is  purely  a  function  ofthe  geometric  boundary  conditions  (captured  by  the  active  constraint  vector, 
{ua},  and  the  corresponding  orthogonal  sorting  matrix,  [ S ])  and  the  two  shape  parameters  ai  and  02.  The  issue  then, 
is  to  deduce  the  active  constraints. 


{u“}|=  f{F  }) 
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Deducing  Active  Constraints 

Whether  a  new  active  constraint  must  be  added  and,  if  so,  which  active  constraint  must  be  added  is  based  on  the 
impingement  of  the  road  surface  on  the  tire  surface.  It  is  tempting  to  simply  compare  the  relative  positions  of  the 
undeformed  shape  of  the  tire  and  the  road  surface  and  sort  the  points  by  the  magnitude  of  impingement.  However; 
one  point  alongthe  road  surface  may  impinge  farther  into  the  undeformed  tire’s  periphery  than  another,  yet  the  former 
point  may  not  be  an  active  constraint  while  the  later  point  is .  This  is  demonstrated  by  the  example  in  Figure  2.  The 
undeformed  tire  (shown  as  a  dashed  line  in  Figure  2a)  is  impinged  uponby  three  distinct  points  alongthe  road  surface 
(shown  as  a  solid  line).  The  points  are  labeled  in  Figure  2b  according  to  their  magnitude  of  impingement,  frommost 
to  least,  as  A,  B,  and  C.  Although  the  impingement  by  pointB  is  greater  than  that  ofpoint  C,  it  should  notbe  identified 
as  an  active  constraint  lest  the  unidirectional  boundary  condition  be  violated  and  the  road  is  presumed  to  “puff’  the 
tire  down.  In  this  way  the  process  cannot  be  based  on  a  Lagrangian  formulation,  since  that  assumes  some  a  priori 
knowledge  of  which  degrees  of  freedom  are  active,  and  thus  motivates  the  development  of  a  recursive  process  for 
adding  active  constraints. 


Figure  2.  Deducing  Active  Constraints 

Although  the  recursive  process  to  identify  the  active  constraints  at  the  tire -road  interface  developed  in  this  study  is 
applicable  to  a  wide  range  oftire  models,  including  FEmodels  and  ring  models,  it  is  described  herein  for  the  tie 
model  defined  by  Equation  (20),  and  shown  schematically  in  Figure  1.  The  first  trial  constraint  mode  is  the 
undeformed  shape  of  the  tire  (a  circle,  partly  shown  as  a  dashed  line  in  Figure  3a),  located  at  a  prescribed  vertical 
positionabovetheroadsurface.  The  recursive  process  continues  through  steps  athrougheas  follows: 

a.  Identify  the  set  of  degrees  of  freedom  where  the  road  surface  impinges  on  the  current  trial  constraint 
mode. 

b.  Sort  the  degrees  of  freedom  with  respect  to  the  magnitude  of  impingement. 

c.  End  the  iteration  process  if  no  significant  impingement  exists. 

d.  Augment  the  current  set  of  active  degrees  of  freedomwith  the  degree  of  freedom  that  has  the 
greatest  iirpingement  (and  those  that  impinge  within  some  very  small  tolerance  of  that  greatest 
inpingement,  if  applicable). 

e.  Develop  a  new  trial  constraint  mode  based  on  this  augmented  set  of  active  constraints. 

For  clarity,  the  set  of  active  constraints  at  the  tire -road  interface  is  developed  through  a  typical  example.  Consider  a 
tire  on  a  flat  plate  as  shown  in  Figure  3  where  the  road,  undeformed  shape  of  the  tire,  and  the  deformed  shape  (the 
trial  constraint  mode)  are  shown  as  a  solid,  dashed,  and  dotted  line  respectively.  The  undeformed  tire  is  the  first  trial 
constraint  mode,  shown  in  Figure  3a.  The  tire  segment  with  the  largest  inpingement  is  directly  belowthe  tire  center 
and  becomes  the  first  active  degree  of  freedom  Figure  3b  shows  the  second  trial  constraint  mode  and  the  first  active 
degree  of  freedomis  shown  as  a  single  small  circle.  Figure  3c  shows  that  two  degrees  of  freedomhave  the  same 
impingement  (due  to  symmetry)  andboth  are  added  to  the  set  of  active  constraints  and  a  new  trial  constraint  mode  is 
calculated.  This  process  ends  with  Figure  3d  when  no  significant  impingement  remains . 

\  \  / 
c  d 

Figure  3.  Recursive  process  to  identify  active  constraints 

Nonlinearity  from  Unidirectional  Geometric  Boundary  Condition 

It  is  clear  that  the  stiffness  equation,  developed  as  Equation  (8)  for  the  continuous  case  and  Equation  (11)  for  the 
discrete  case,  are  linear  and  time  invariant  as  are  the  physical  properties  shown  in  Table  1.  However  real  tires  exhibit 
a  nonlinear  vertical  force-deflection  behavior  that  must  be  captured.  The  unidirectional  geometric  boundaiy 
conditions  provide  the  mechanism  by  which  this  non-linear  stiffness  relationship  is  manifest.  As  the  tire  deflection 
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increases,  more  tire  segments  come  into  contact  with  the  road  surface  and  become  active  constraints  as  shown  in 
Figure  4.  As  the  tire  first  touches  the  road  surface,  only  two  segments  are  active  and  produce  forces,  shown  as  arrows 
in  Figure  4a.  As  the  deflection  increases,  the  magnitude  of  the  forces  at  the  existing  active  constraints  increases  and 
additional  degrees  of  freedombecome  active,  as  shown  by  the  four  arrows  in  Figure  4b.  This  addition  of  new 
constraints  and  the  increase  force  magnitude  continues  with  increased  deflection  (shown  in  Figure  4c). 


a  b  c 


Figure  4.  Additional  radial  forces  added  as  deflection  increases 


EXPERIMENTAL  VALIDATION 

The  experimental  data  used  to  validate  this  work  is  provided  by  Professor  Schalk  Els  of  the  University  of  Pretoria, 
South  Africa.  The  tire  used  is  a  Continental Conti-Trac  AT238/85  R16.  Two  quasi-static  tests  (a  flat  plate  test  anda 
cleat  test)  are  performed  that  produce  a  resultant  spindle  force,  F,  with  respect  to  the  tire  deflection,  e ,  as  shown  in 
Figure  5a  for  the  flat  plate  test  andFigure  5b  for  the  cleat  test.  The  side  lengthofthe  square  cleat  is  19  mm  In  both 
cases  an  obstruction  (plate  or  cleat)  is  gradually  pressed  toward  the  center  of  a  tire  that  is  rigidly  fixed  at  the  spindle. 
The  model  stiffness,  ko ,  and  model  shape  parameters,  ai  and  a?,  are  identified  using  these  two  quasi-static  tests.  The 
flat  plate  test  results  are  compared  to  the  results  fromRagheb,  El-Gindy,  and  Kishawy  [7].  The  results  ofthese  two 
tests  are  shown  in  Figure  6  and  Figure  7;  in  both  figures  the  experimental  results  are  shown  as  a  solid  line  while  the 
simulation  results  are  shown  as  a  dashed  line. 


a.  Flat  Plate 


b.  Cleat 


Figure  5.  Schematics  for  Test  Procedures 


Validation  via  Flat  Plate  Testing 

It  is  clear  from  Figure  6  that  the  nonlinear  stiffening  characteristic  of  the  tire  is  capturedby  the  constraint  mode  tie 
model  via  the  unidirectional  geometric  boundary  condition  (see  also  Figure  4).  The  model  parameters  were  optimized 
to  minimize  the  maximum  error,  in  this  case  4%,  throughout  a  reasonably  large  range  of  travel.  For  comparison 
purposes,  Ragheb,  El-Gindy,  and  Kishawy  [7]  used  Mooney-Rivlin  material  properties  to  create  detailed  FE  models 
of  the  tread  and  sidewall  for  a  tire  similar  to  the  Michelin  Off-road  12.00R20  XMLTL  149J  (a  somewhat  larger  tire 
than  that  used  in  the  present  study).  The  FE  model  produced  quasi-static  flat -plate  test  results  that  were  within  10- 
15%  of  the  experimental  results  (depending  on  tire  pressure).  Also  ofnoteis  that  the  nonlinearity  in  their  simulation 
was  negligible  whereas  the  current  study  produces  results  that  are  clearly,  thoughnot  dramatically,  nonlinear. 
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Deflection  (cm) 


Figure  6:  Flat-plate  experimental  and  simulated  force,  F,  with  respect  to  the  deflection,  e 

In  this  example,  the  shear  and  radial  stiffnesses  are  orders  ofmagnitude  smaller  than  the  bending  stiffness.  To  simplify 
future  comparisons  to  these  results,  the  model  shape  and  stiffness  parameters  and  tire  stiffh  ess  parameters  provided  in 
Table  2  may  be  used. 


Table  2.  Simplified  tire  stiffness  parameters  and  modeling  parameters 


Modeling 

Parameters  (N=72) 

Tire  Stiffness  Parameters 

fcn  =  7.053  x  106 

Yi  =3.15 

a ,  =  -2/3 

O 

II 

£ 

cc  2  — 1/6 

7s  =  0 

Validation  via  Cleat  Testing 

The  cleat  test  simulation  results  shown  in  Figure  7  show  very  good  agreement  throughout  the  range  of  deflection;  the 
maximum  error  is  7%.  The  nonlinear  stiffness  property,  in  the  form  of  a  discontinuity,  is  clearly  captured  in  this 
modeling  technique.  The  discontinuity  in  stiffness  occurs  when  the  experimental  deflection  exceeds  approximately 
5cm,  dividing  the  results  into  two  regions.  In  the  lower  region  the  tire  is  suspended  by  the  cleat  and  the  contact  aiea 
is  constrained  to  the  cleat  surface  only.  In  this  lo wer  region  the  stiffness  is  quite  linear  in  both  the  experimental  results 
and  the  simulated  results.  This  is  ejected  because  there  are  no  additional  unidirectional  boundary  constraints  added 
as  the  deflection  increases  in  this  lower  region.  There  is  a  dramatic  change  in  the  stiffness  as  deflection  increases  to 
the  point  where  the  tire  tread  touches  the  surface  below  the  cleat  as  shown  in  the  lower  inset  picture  in  Figure  7.  Note 
that  this  discontinuity  in  stiffness  is  captured  by  the  constraint  mode  tire  model  even  though  a  simple  linear  model  is 
used  to  predict  the  radial  deflection  of  the  tire  surface. 
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Figure  7:  Cleat  experimental  and  simulated  force,  F,  with  respect  to  the  deflection,  e 


DISCUSSION 

It  is  envisioned  that  the  constraint  mode  tire  model  developed  in  this  work  canserveas  a  morphological  filter.  That 
is,  this  constraint  mode  tire  model  can  be  used  to  pre -filter  the  surface  once,  providing  the  required  bridging  and 
enveloping  properties  ofthe  tire.  A  simpler  tire  model  (perhaps  a  simple  linear  point -follower)  could  be  usedin  the 
iterative  design  process  to  provide  fast  yet  reliable  spindle  force  predictions  for  vehicle  dynamic  simulation  and 
reliability  evaluation.  Specifically,  an  effective  roadprofile  could  be  estimated  given  an  actualprofile  by  simulating 
the  constraint  mode  tire  model  traveling  over  an  actual  profile  given  a  constant  spindle  load.  It  is  hoped  that  a  small 
sacrifice  in  accuracy  will  result  in  an  order  of  magnitude  increase  in  computational  speed.  Preliminary  results  shown 
in  Figure  8  where  the  actualprofile  is  shown  as  a  solid  grey  line  and  the  effective  road  profile  is  shown  as  a  dashed 
black  line.  In  this  example  both  the  bridging  and  enveo loping  proerties  are  evident  (including  the  bridging  of  deep 
narrow  cracks)  which  is  encouraging  for  further  study. 


CONCLUSIONS 

Themain  contribution  of  this  work  is  the  development  ofa  computationally  efficient,  linear,  planar,  quas  i -static  tie 
model  that  accurately  predicts  a  tire’s  lower-frequency,  reasonably  large  amplitude,  nonlinear  stiffness  relationship. 
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One  major  development  is  that  themodelis  parameterized  by  a  single  stiffness  parameter  and  two  shape  parameters. 
The  implication  of  this  development  is  that  the  tire’s  deformed  shape  is  completely  independent  of  the  overall  tiie 
stiffness  and  the  forces  acting  onthetire.  Rather,  the  shape  is  purely  a  functionofthe  geometric  boundary  conditions 
and  the  two  shape  parameters  ai  and  0L2.  This  simple  model  produces  nonlinear  stiffiiess  when  simulating  a  flat  plate 
test  and  a  discontinuous  stiffiiess  when  simulating  a  cleat  test,  in  both  cases  producing  reasonable  force  predictions. 
The  nonlinear  behavior  stems  from  the  unidirectional  geometric  boundary  conditions  that  are  preserved  during  the 
recursive  method  developed  to  deduce  the  active  constraints.  A  simple  example  demonstrates  that  simpler  algorithms 
that  do  notusethis  new  recursive  method  can  erroneously  define  non-active  constraints  tobeactive.  This  may  produce 
the  effect  oftheroadsurface  “pulling”  the  tire  down.  This  model  strikes  a  balance  between  heuristic  tire  models  (such 
as  a  linearpoint-follower)  thatlackthe  fidelity  to  make  accurate  chassis  load  predictions  and  computationally  intensive 
models  that  cannot  provide  timely  predictions.  It  is  hoped  that  this  tire  model  is  used  as  an  integral  part  of 
computationally  efficient  and  accurate  vehicle  dynamic  simulations  that  are  critical  throughout  the  iterative  design 
process. 
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